Parabolic cylinder functions implemented in Matlab 

E. Cojocaru 
Department of Theoretical Physics, Horia Hulubei 
National Institute of Physics and Nuclear Engineering, 
Magurele- Bucharest P.O.Box MG-6, 077125 Romani^ 

Routines for computation of Weber's parabolic cylinder functions and their deriva- 
fives are implemented in Matlab for both moderate and great values of the argument. 
Standard, real solutions are considered. Tables of values are included. 
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I. INTRODUCTION 

The parabolic cylinder functions were introduced by Weber [T] in 1869. Standard solu- 
tions to Weber's equation were given by Miller [2] in 1952. These relations are also provided 
by Abramowitz and Stegun [3]. There are two standard forms of the Weber's equation, 

g-(|x 2 + a), = 0, (1) 

g + {\x* - a)y = 0. (2) 

Equation ([2J is obtained from ([!]) with changes a by — ia and x by xe* 71 "/ 4 . Thus, if y(a, x) 
is a solution of (fl|, then ^ has solutions: y{— ia, xe J7r//4 ), y(— ia, — xe l7f//4 ), y(ia, — xe - "^ 4 ), 
and y(ia, xe~ l7T / 4 ). In the following we consider only real solutions of real equations. 



II. SOLUTIONS OF EQUATION ([!]) 

A. Standard solutions 

There are two standard solutions of Eq. ([I]), U(a, x) and V(a, x), both of them expressed 
in terms of Whittaker's function D n _i, 

a 2 

U(a,x) = D_ a _i, (3) 

V(a,x) = ^T(l + a)[sm(7ra)U(a,x) + U(a, -x)]. (4) 
In a more symmetrical notation, these solutions are 

U(a,x) = D_ a _i_ = Y lC osp -Y 2 smp, (5) 

V(a,x) = * (Fisin/3 + r 2 cos/j), (6) 
L [ 2 a) 



where 



/3 = tt(| + |), (7) 

y _ vin\ - 1 ) y _ yjni - 1 ) , 8 . 

^^2 2^4 v /7T2 2 4 

, x 

y 1 = l + a- + (a 2 + -)— + (a 3 + -)- + (a 4 + 11a 2 + —)— + ..., (9) 

V2=x + a^ + (a 2 + ^ + (a 3 + + (a 4 + 17a 2 + ^)|- + . . . , (10) 



in which the coefficients A n of ^ obey the recurrence relation 

A n+2 = aA n + \n(n - l)A n _ 2 - 

Similarly to Eq. (|4]), there is relation 

U (a, x) = — * — - — -\V(a, -x) - sin(7fa)V{a, x)]. 
1 {a + |J cos z (7ra) 

At x = 0, 



U(a,0) 



U'(a,0) 



V{a,0) 



2 f+ir(f + |)' 

7T 

1 



2f-ir(§ + 1 

21 +3 sin7r 



4 2- 
3 a \ 



r( 4 2 

2f+f SHITT^ - ^ 



r(l-f) 

B. Recurrence relations for U(a,x) and ^(a,x) 

Standard solutions {/(a, x) and V(a, x) obey the recurrence relations 

xU(a, x) — U(a — 1, x) + (a + \)U(a + 1, x) = 0, 
U'(a, x) — \x\J (a, x) + U(a — 1, x) = 0, 

xV(a,x) - V{a + l,x) + (a - \)V{a -l,x) = 0, 
V'(a, x) — ^xV(a, x) — (a — \)V(a — 1, x) = 0. 

C. Relations at large values of argument x 

At large values of argument x, when x 3> \a\, there are relations 

| (a + |)(a + f ) (a + |)(a + |)(a + §)(a + |) 



[/(a, x) ~ x a 2 e 



1 + 



2x 2 2 • 4x 4 

(a-f)(a-§) , (o-|)(o-|)(o-|)(a-|; 



2x 2 2 ■ 4x 4 
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D. Analytic relations at specific values of parameter a 

At half of odd parameter a, the standard solutions obey relations given in Table [Tj 
TABLE I: Analytic relations for U(a, x) and V(a, x) at half of odd parameter a 



U{- 


-0.5, x) 


X 2 

= e 4 






t/(o.5,x) = A 


1 ey 3p" 

/ 4 

V 


U(- 


-1.5, x) 


= xe 4 






y(i.5,x) = y 


/ -xe 4 

v 77 


U(- 


-2.5, x) 


= (x 2 - 


l)e~T 




V(2.5,x) = A 


Jl(x 2 + l)e^ 


U{- 


-3.5, x) 


= (x 3 — 


3x)e 4 




y(3.5,x) = y 


/f (x 3 + 3x)e^ 


U(- 


-4.5, x) 


= (x 4 - 


6x 2 + 3)e _ 


.2 
4 


V(4.5,x) = y 


/J (x 4 + 6x 2 + 3)eT 



There is also relation, 



[7(0.5, x) = «/§e^erfc( 



(19) 



2 ~"~ V V2 ; 

where erfc(;s:) is the complementary error function [3]. Using recurrence relations ( Jl5| ) and 
taking into account that 

d erfc(4?) /— j2 

(20) 



one obtains 



Further we obtain 



V2' 



dx 



2 

z e 2 , 



7'(0.5,x) = |«7(0.5, x) - e~T 



[7(1.5, x) = -x7(0.5,x) + e" 
7(2.5, x) = ^(x 2 + 1)7(0.5, x) - \xe 



(21) 



(22) 



7'(1.5,x) = -(y + l)7(0.5,x) + |e-T, 
7'(2.5,x) = |x(x 2 + 5)7(0.5, x) - ( y + l)e - ^. 



(23) 



At integer values of parameter a, the standard solutions 7(a, x) and V(a,x) obey relations 
given in Table [Lt] They are expressed in terms of modified Bessel functions I v (z) and K„(z), 
with 



I-u - Iy = l sin(-Ku)K v , 



I-v + h = cos(iris)3 u . 



(24) 
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TABLE II: Analytic relations for U(a,x) and V(a,x) at integer values of parameter a 



U(l,x) = 27T-2(f )§(-#! +Kz) 

4 4 

U(2,x) = 1^-^(1)1(2^1 - 3Kt +Ks) 

4 4 4 


V(0,x) = 

F(l,x) = Ifipi -2fs) 

4 4 

F(2,x) = if i(23Ti -3J 3 + J 5 ) 

z z 4 4 4 


U{-l,x) =7r-3(f)i(jri +ifs) 

4 4 

C/(-2,x) =^-2(1)1(2^1 +3K 3 -fTs) 

4 4 4 y 


y(-i,x) = f ipi + dt 3 ) 

4 4 

F(-2,x) = |f l(2Ji +3^3 -3s) 

13 Z A A A' 



The argument of all Bessel functions in Table II is \. Further relations can be obtained by 
using the recurrence relations for the modified Bessel functions, 



ZKI(Z) : 



-zK v Az) - vKJz) 



(25) 



2v 

Iu+l{z) = L-l{z) - —I V (Z), 

zl' v (z) = zl v -i(z) - vl v [z), 
with K- V (z) = K u (z) and d- v (z) = 3 v (z). 



(26) 



III. SOLUTIONS OF EQUATION [5] 



A. Standard solution 



The standard solution W(a,x) of Eq. ([2| is 



W(a,±x) = 2-*(\l- r y lT 

^3 




2G, 
~G~i 



where 



G 1 = m + i 



yi 



V2 



1 + a— + (a 2 - -)— (a 



4715 
3 



.r 



2' 4! 

3,:r 5 



x + a— + (a -)— + {a J 



G 3 = |r(f + |)|, 

y)^ + (a 4 -Ha 2 + -) 

X3& , . a r> 63 , x 



+ . . 



+ 



in which the coefficients A n of obey the recurrence relation 



.4 



n+2 



\n{n - l)A»-2- 



(27) 

(28) 
(29) 

(30) 
(31) 
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Relations for Gamma function of complex argument are given in Appendix [Aj At x = 0, 




W(o,0) = 2-*y G , . 
W'(a,0) = -2-^f^. (32) 

V Lri 



B. Relations at large values of argument x 



At large values of the argument x, when x 3> |a|, there are relations 

l~2k 

W(a,x) = a/ — [si(a,x) cos7 — s 2 (a,x) sin 7], 



W(a, —x) 



y^[si(a,x) sin 7 + s 2 (a,x) cos 7], (33) 



where 



with 



with 



A; = Vl + e 2 ™ - e™, AT 1 = Vl + e 2 ™ + e™, (34) 
7= X _a 4 + 2' ^ 



= argT(za + |), (36) 

/ x 1 . ^2 «4 ^6 U 8 . . 

SlM ~ 1+ l!tf-2!2V-3!2¥ + 4!2¥ + "-' (3?) 

S2(a ' X) ~ ~ T!2^ " 2!2^ + i¥ + i¥ " " " ' (38) 
T(m + ia + |) 

u m + iVm = — ^r. — ' Y\ — > m = 2,4,6... (39) 



C. Analytic relations at a = 

At a = there are relations 

M/(0,±x) = 2-*ViF£[J_j(^) =F Jj(^)], (40) 

oW(0, ±x) 9 , — rT .x 2 . .x 2 . . 

d ; = -2"^v^[J f ( T ) ± J_ f ( T )], (41) 

where J v {z) is the Bessel function of the first kind. 
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IV. IMPLEMENTATION OF PARABOLIC CYLINDER FUNCTIONS IN 

MATLAB 

Routines implemented in Matlab [4J for computation of parabolic cylinder functions are 



shortly described in Table III For moderate values of argument x and parameter a, standard 
parabolic cylinder functions U(a, x) and V(a, x) are computed with routines "pu" and "pv", 
respectively, whereas W(a,x) is computed with routine "pw". Differentiation with respect 
to argument x is computed with routines "dpu", "dpv", and "dpw". For large values of 
argument x, when |x| 3> \a\, functions U(a,x), V(a, x), and W(a,x) are computed with 
routines "pulx" , "pvlx" , and "pwlx" , and their derivatives with routines "dpulx" , "dpvlx" , 
and "dpwlx", respectively. Routine "cgamma" computes the gamma function of complex 
argument; using the function code kf, it computes either the logarithm of gamma function 
(when kf = 0) or gamma function (when kf = 1). Values of parabolic cylinder functions 



obtained by using these routines are shown in Tables IV IX 



APPENDIX A: RELATIONS FOR GAMMA FUNCTION OF COMPLEX 

ARGUMENT 

If \z\ 3> 1 and | &rgz\ < ir — e with e > 0, there is relation 

°° B 1 

lnr(z) ~ (z - |) In, - z + |ln 27 r + ^ 2n{2 ^_ x) z2n -v ( A1 ) 
where B 2n are the Bernoulli's numbers, 

B* = (-ir i2 (Mf:^> *=1,2,... (A2) 

^ ' n=l 

Specific values, 



B 2 


= 1/6 


Bu 


= -691/2730 


B 4 


= -1/30 


B\4 


= 7/6 


B 6 


= 1/42 


Bw 


= -3617/510 


B 8 


= -1/30 


Big 


= 43867/798 


Bio 


= 5/66 


B20 


= -174611/330 



(A3) 
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Other useful relations, 



r(z + n) = z(z + 1) ■ ■ ■ (z + n - l)V(z) 



mn-z) = —^- v (A4) 

z sin (7TZ) 



[1] H. F. Weber "Ueber die Integration der partiellen Differential-gleichung: d 2 u/dx 2 + d 2 u/dy 2 + 

k 2 u = 0," Math. Ann. 1, 1-36 (1869). 
[2] J. C. P. Miller "On the choice of standard solutions to Weber's equation," Proc. Cambridge 

Philos. Soc. 48, 428-435 (1952). 
[3] M. Abramowitz and I. Stegun Handbook of Mathematical Functions (New York, 1964). 



[4] E. Cojocaru, Matlab free available computer code Parabolic cylinder functions 
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TABLE III: Routines for parabolic cylinder functions 



Name of routine 



Routine call 



What the routine computes 



cgamma 



pu 



dpu 



pv 



dpv 



pw 



dpw 



pulx 



dpulx 



pvlx 



dpvlx 



pwlx 



dpwlx 



[gr, gi] =cgamma(a;, y, kf) 



u=pu(a, x) 
du=dpu(a, x) 
v=pv(a, x) 
dv=dpv(a, x) 
w=pw(a, x) 
dw=dpw(a, x) 
it=pulx(a, x) 
du=dpulx(a, x) 



w=pvlx(a, x) 
ofc=dpvlx(a, x) 



w=pwlx(a, x) 
dw=dpw1x(a, x) 



r(z) with complex argument z (when kf = 1) or lnT(z), 
(when kf = 0); x and y are the real and imaginary parts of z 
gr and gi are the real and imaginary parts of lnT(z) or T{z) 



[Eqs. (Al 



A4)]. 



Parabolic cylinder function U(a,x) for moderate values of 
parameter a and argument x [Eqs. (3-12)]. 



Derivative with respect to x of parabolic cylinder function 
U(a, x) for moderate values of parameter a and argument x 
Parabolic cylinder function V(a, x) for moderate values of 
parameter a and argument x [Eqs. (3-12)]. 



Derivative with respect to x of parabolic cylinder function 
V(a, x) for moderate values of parameter a and argument x 
Parabolic cylinder function W(a, x) for moderate values of 



3D] 



parameter a and argument x [Eqs. (27- 
Derivative with respect to x of parabolic cylinder function 
W(a, x) for moderate values of parameter a and argument x 
Parabolic cylinder function U(a,x) for large values of parameter 



x (\x\ 3> \a\) and moderate values of parameter a [Eq. (17)]. 
Derivative with respect to x of parabolic cylinder function 
U(a,x) for large values of parameter x (\x\ 3> \a\) and moderate 
values of parameter a. 

Parabolic cylinder function V(a, x) for large values of parameter 



x {\x\ 3> \a\) and moderate values of parameter a [Eq. (18)]. 
Derivative with respect to x of parabolic cylinder function 
V(a,x) for large values of parameter x (\x\ 3> \a\) and moderate 
values of parameter a. 

Parabolic cylinder function W(a, x) for large values of parameter 



39)] 



x {\x\ S> |a|) and moderate values of parameter a [Eqs. (33 
Derivative with respect to x of parabolic cylinder function 
W(a,x) for large values of parameter x (\x\ 3> \a\) and moderate 



TABLE IV: Values of U(a, x) with a = -5, -3.5, -1, 1, 3.5, 5 and x = 0, 1, 3, 5 



x\q 


-5.0 


-3.5 


-1.0 


0.0 


3 052183664350372 


-0 000000000000000 


581368317019118 


1.0 


0.579926011661105 


-1.557601566142810 


0.842203244069839 


3.0 


3.202129097812791 


1.897186042113549 


0.184881790005045 


5.0 


1.879976816310843 


0.212349954984646 


0.004337473181400 


x\a 


1.0 


3.5 


5.0 


0.0 


1.162736634038237 


0.333333333333333 


0.103354367470066 


1.0 


0.378262434740955 


0.048971230815929 


0.010659966828235 


3.0 


0.017224293634316 


0.000610423938072 


0.000070950238455 


5.0 


0.000161381143270 


0.000002208878109 


0.000000155227075 



TABLE V: Values of U{a, -x) with a = -5, -3.5, -1, 1, 3.5, 5 and x = 0, 1, 3, 5 



x\a 


-5.0 


-3.5 


-1.0 


0.0 


3.052183664350372 


-0.000000000000000 


0.581368317019118 


1.0 


-4.332232266251285 


1.557601566142810 


-0.195001018223362 


3.0 


3.802753160685226 


-1.897186042113549 


-1.767855400724101 


5.0 


-9.615606269532364 


-0.212349954984649 


-35.754085404247576 


x\a 


1.0 


3.5 


5.0 


0.0 


1.16273663404 


0.33333333333 


0.10335436747 


1.0 


3.27078479478 


2.19468750736 


0.97838806074 


3.0 


45.73101176423 


142.69397188181 


125.30190015651 


5.0 


3259.12460949910 


30297.53050402874 


45998.28922772748 



TABLE VI: Values of V(a, x) with a = -5, -3.5, -1, 1, 3.5, 5 and x = 0, 1, 3 





-5.0 


-3.5 


-1.0 


0.0 


-0 058311457540778 


265961520267622 


-0 656003897333753 


1.0 


0.082766571619165 


-0.076762147625440 


0.220035086525655 


3.0 


-0.072650962016911 


0.097154672861824 


1.994811204614366 


5.0 


0.183704546768818 


1.173350875864019 


40.344165108706711 


x\a 


1.0 


3.5 


5.0 


0.0 


0.3280019487 





1.7220102305 


1.0 


0.9226713556 


4.0980162226 


16.3011422859 


3.0 


12.9004802412 


272.5242458690 


2087.6829809173 


5.0 


919.3820780818 


57864.0209141053 


766387.7838412275 



TABLE VII: Values of V(a, -x) with a = -5, -3.5, -1, 1, 3.5, 5 and x = 0, 1, 



x\a 


-5.0 


-3.5 


-1.0 


0.0 


-0.058311457540778 


0.265961520267622 


-0.656003897333753 


1.0 


-0.011079389291262 


-0.076762147625440 


-0.950324595068664 


3.0 


-0.061176139925034 


0.097154672861824 


-0.208616760217021 


5.0 


-0.035916642101972 


1.173350875864019 


-0.004894314375732 


x\a 


1.0 


3.5 


5.0 


0.0 


0.32800194867 





1.72201023050 


1.0 


0.10670586276 


-4.09801622261 


0.17760809131 


3.0 


0.00485888353 


-272.52424586904 


0.00118211779 


5.0 


0.00004552478 


-57864.02091410524 


0.00000258678 



TABLE VIII: Values of W(a, x) with a = -5, -3, -1, 1, 3, 5 and x = 0, 1, 3, 5 



x\(X 


-5.0 


-3.0 


-1.0 


0.0 


473478576486605 


539330386270653 


0.731481090245431 


1.0 


-0.657520526362908 


-0.611126375982879 


-0.184115556183355 


3.0 


-0.062604004232077 


0.636305300554784 


-0.053352644054153 


5.0 


0.089361847055232 


0.437066960213013 


-0.570254174032845 


x\a 


1.0 


3.0 


5.0 


0.0 


0.731481090245431 


0.539330386270653 


0.473478576486605 


1.0 


0.315937643962764 


0.101682226485666 


0.052572013487910 


3.0 


0.016773032899024 


0.009166528652640 


0.001223742332881 


5.0 


0.022807516888135 


-0.003844865237560 


0.000115773464320 



TABLE IX: Values of W(a, -x) with a = -5, -3, -1, 1, 3, 5 and x = 0, 1, 3, 5 



x\a 


-5.0 


-3.0 


-1.0 


0.0 


0.473478576486605 


0.539330386270653 


0.731481090245431 


1.0 


0.070610950611453 


0.428801301530536 


0.950916920458344 


3.0 


0.606270877302830 


0.177268761402591 


-0.757374330077355 


5.0 


0.538608396875686 


-0.370945283780393 


0.180907184885679 


x\a 


1.0 


3.0 


5.0 


0.0 


0.731481090245 


0.539330386271 


0.473478576487 


1.0 


1.903689596383 


3.001251077335 


4.378212848013 


3.0 


6.183176599808 


57.210355295947 


253.398744868662 


5.0 


-4.359927574948 


66.590129609337 


2852.835947866653 



